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Searches for lepton flavour and lepton number violation in kaon decays by the NA48/2 and 
NA62 experiments at CERN are presented. A new measurement of the helicity suppressed 
ratio of charged kaon leptonic decay rates Rk = T(K e 2) /T^K^) to sub-percent relative pre- 
cision is discussed. An improved upper limit on the lepton number violating — > 7r T fi ft 
decay rate is also presented. 



1 Introduction 

In the Standard Model (SM) the decays of pseudoscalar mesons to light leptons are helicity 
suppressed. In particular, the SM width of P ± — >• i^v decays with P = it, K, D, B (denoted P12 
in the following) is 



where Gp is the Fermi constant, Mp and are meson and lepton masses, fp is the decay 
constant, and V qq! is the corresponding Cabibbo-Kobayashi-Maskawa matrix element. Although 
the SM predictions for the P12 decay rates are limited by hadronic uncertainties, their specific 
ratios do not depend on fp and can be computed very precisely. In particular, the SM pre- 
diction for the ratio Rk = T{K e 2) /T{K^2) of kaon leptonic decay widths inclusive of internal 
bremsstrahlung (IB) radiation is^ 

RT = (^) 2 {^3) 2 C 1 + = ( 2 - 477 ± °- 001 ) X 1(r5 ' ( 2 ) 

where 5Rqed = (—3.79 ± 0.04)% is an electromagnetic correction due to the IB and structure- 
dependent effects. 

Within certain two Higgs doublet models (2HDM of type II) , including the minimal super- 
symmetric model (MSSM), Rk is sensitive to lepton flavour vio latin g (LFV) effects appearing at 
the one- loop level via the charged Higgs boson (H ^ ) exchangers, representing a unique probe 
into mixing in the right-handed slepton sector^. The dominant contribution due to the LFV 



coupling of the is 



(3) 

where tan/3 is the ratio of the two Higgs vacuum expectation values, and |A^| is the mixing 
parameter between the superpartners of the right-handed leptons, which can reach ~ 10 -3 . 
This can enhance Rk by 0(1%) without contradicting any experimental constraints known at 
present, including upper bounds on the LFV decays r — > eX with X = rj, 7, fifi. On the other 
hand, Rk is sensitive to the neutrino mixing parameters within the SM extension involving a 
fourth generation^. 

The first measurements of Rk were performed in the 1970s EES while the current PDG 
world average^ is based on a more precise recent result^ R K = (2.493 ± 0.031) x 10" 5 . A 
new measurement of Rk based on a part of a dedicated data sample collected by the NA62 
experiment at CERN in 2007 is reported here: the analyzed K e 2 sample is ~ 4 times larger than 
the total world sample, allowing the first measurement of Rk with a relative precision below 
1%. 

The decay K — > 7r =F ^ ± /i = ' = violating lepton number by two units can proceed via a neutrino 
exchange if the neutrino is a Majorana particle, consequently the experimental limits on this 
decay provide constraints on the effective Majorana neutrino mass (%)^. This decay has 
also been studied in the context of supersymmetric models with i?-parity violation^]. The best 
previous upper limit on the decay rate was based on a special data set collected by the BNL 
E865 experiment in 1997031 The sample of 717/// triggers collected by the NA48/2 experiment at 
CERN during the 2003-04 data taking is about 8 times larger than the E865 one, which allows 
improving the upper limit significantly. 
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2 Beam and detector 

The NA48/2 and NA62 (phase I) experiments at CERN took data in 2003-04 and 2007-08, 
respectively, using the same kaon beamline and experimental setup^l The trigg er lo gic was 
optimized to detect direct CP violating charge asymmetries in decays in 2003-04^21, and for 
the K&jKpfi ratio measurement in 2007-08. The beam line is capable of delivering simultaneous 
unseparated K + and K~ beams derived from the 400 GeV/c primary proton beam extracted 
from the CERN SPS. Central values of kaon momentum of 60 GeV/c (both K + and K~ beams) 
and 74 GeV/c (K + beam only), with a narrow momentum band, were used for collection of the 
main data samples by the NA48/2 and NA62 experiments, correspondingly. 

The fiducial decay region is contained in a 114 m long cylindrical vacuum tank. With ~ 10 12 
primary protons incident on the target per SPS pulse of 4.8 s duration, the typical secondary 
beam flux at the entrance to the decay volume is 10 7 to 10 8 particles per pulse, of which about 
5% are kaons, while pions constitute the dominant component. The transverse size of the beams 
within the decay volume is below 1 cm (rms), and their angular divergence is negligible. 

Among the subdetectors located downstream the decay volume, a magnetic spectrometer, a 
plastic scintillator hodoscope (HOD), a liquid krypton electromagnetic calorimeter (LKr) and a 
muon veto counter (MUV) are principal for the present measurements. The spectrometer, used 
to detect charged products of kaon decays, is composed of four drift chambers (DCHs) and a 
dipole magnet. The HOD producing fast trigger signals consists of two planes of strip-shaped 
counters. The LKr, used for particle identification and as a veto, is an almost homogeneous 
ionization chamber, 27Xo deep, segmented transversally into 13,248 cells (2x2 cm 2 each), and 
with no longitudinal segmentation. The MUV is composed of three planes of plastic scintillator 
strips read out by photomultipliers at both ends. A beam pipe traversing the centres of the 



detectors allows undecayed beam particles and muons from decays of beam pions to continue 
their path in vacuum. 



3 Search for lepton flavour violation 

The precision measurement of Rk = T(K e2 ) /TiK^) is based on the NA62 2007 data sample. 
The measurement method is based on counting the numbers of reconstructed K e2 and K^ 2 can- 
didates collected concurrently. Consequently the result does not rely on kaon flux measurement, 
and several systematic effects (e.g. due to reconstruction and trigger efficiencies, time-dependent 
effects) cancel to first order. To take into account the significant dependence of signal accep- 
tance and background level on lepton momentum, the measurement is performed independently 
in bins of this observable: 10 bins covering a lepton momentum range of (13; 65) GeV/c are 
used. The ratio Rk in each bin is computed as 

1_ N(K e2 ) - N B (K e2 ) A{K^ 2 ) x e(ggg) 1 
K D ' N(K^ 2 ) - iV B (i^ 2 ) ' A(K e2 ) ' f e x e(K e2 ) ' / LKr ' 1 ' 

where N(K^ 2 ) are the numbers of selected Kg 2 candidates (£ = e,fi), N-q(K^ 2 ) are numbers of 
background events, A(K^ 2 ) / A(K e2 ) is the geometric acceptance correction, are efficiencies of 
e//j, identification, e(K^ 2 ) are trigger efficiencies, /LKr is the global efficiency of the LKr readout, 
and D = 150 is the downscaling factor of the K^ 2 trigger. 

A detailed Monte Carlo (MC) simulation including beam line optics, full detector geome- 
try and material description, stray magnetic fields, local inefficiencies of DCH wires, and time 
variations of the above throughout the running period, is used to evaluate the acceptance cor- 
rection A(K^ 2 )/A(K e2 ) and the geometric parts of the acceptances for background processes 
entering the computation of Nb(K^ 2 ). Simulations are used to a limited extent only: particle 
identification, trigger and readout efficiencies are measured directly. 

Due to the topological similarity of K e2 and K^ 2 decays, a large part of the selection con- 
ditions is common for both decays: (1) exactly one reconstructed positively charged particle 
compatible with originating from a beam K decay; (2) its momentum 13 GeV/c < p < 65 GeV/c 
(the lower limit is due to the 10 GeV LKr energy deposit trigger requirement); (3) extrapolated 
track impact points in subdetectors are within their geometrical acceptances; (4) no LKr energy 
deposition clusters with energy E > 2 GeV not associated to the track, to suppress background 
from other kaon decays; (5) distance between the charged track and the nominal kaon beam axis 
CDA < 3 cm, and decay vertex longitudinal position within the nominal decay volume. 

The following two principal selection criteria are different for the K e2 and K^ 2 decays. Kg 2 
kinematic identification is based on the reconstructed squared missing mass assuming the track 
to be a positron or a muon: M^ iss (£) = (Pk — Pe) 2 , where Pk and Pe (£ = e, /u) are the 
four-momenta of the kaon (average beam momentum assumed) and the lepton (positron or 
muon mass assumed). A selection condition — M 2 < M^ dss (£) < M 2 is applied to select K^ 2 
candidates, where M 2 varies between 0.013 and 0.016 (GeV/c 2 ) 2 and M 2 between 0.010 and 
0.013 (GeV/c 2 ) 2 across the lepton momentum bins, depending on M^ iss (£) resolution. Particle 
identification is based on the ratio E/p of track energy deposit in the LKr calorimeter to its 
momentum measured by the spectrometer. Particles with 0.95 < E/p < 1.1 (E/p < 0.85) are 
identified as positrons (muons). 

Kinematic separation of K e2 from K^ 2 decays is achievable at low lepton momentum only 
(p < 35 GeV/c). At high lepton momentum, the K^ 2 decay with a mis- identified muon (E/p > 
0.95) is the largest background source. The dominant process leading to mis-identification of the 
muon as a positron is 'catastrophic' bremsstrahlung in or in front of the LKr leading to significant 
energy deposit in the LKr. Mis-identification due to accidental LKr clusters associated with 
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Figure 1: Left: Mis-identification probability for muons traversing the lead wall, P^ e b , for (E/p) m \ n = 0.95 as a 
function of momentum: measurement (solid circles with error bars) and simulation (solid line). Right: Correction 
factors /pb = P/ie/Pue f° r the considered values of (B/p) m in , as evaluated with simulation. Dotted lines in both 
plots indicate the estimated systematic uncertainties of the simulation. 



the muon track is negligible, as concluded from a study of the sidebands of track-cluster time 
difference and distance distributions. 

The muon mis-identification probability P^ has been measured as a function of momentum. 
To collect a muon sample free from the typical ~ 10 -4 positron contamination due to fi — >• e 
decays, a 9.2Xq thick lead (Pb) wall covering ~ 20% of the geometric acceptance was installed 
approximately 1.2 m in front of the LKr calorimeter (between the two HOD planes) during 
a period of data taking. The component from positrons which traverse the Pb wall and are 
mis-identified as muons from Ku2 decay with p > 30 GeV/c and E/p > 0.95 is suppressed down 
to a negligible level (~ 10~ 8 ) by energy losses in the Pb. 

However, muon passage through the Pb wall affects the measured P^ via two principal 
effects: 1) ionization energy loss in Pb decreases P^ e and dominates at low momentum; 2) 
bremsstrahlung in Pb increases P^e cind dominates at high momentum. To evaluate the correc- 
tion factor /pb = P^ e /P^e, a dedicated MC simulation based on Geant4 (version 9.2) has been 
developed to describe the propagation of muons downstream from the last DCH, involving all 
electromagnetic processes including muon bremsstrahlung^^. 

The measurements of Pj^ in momentum bins compared with the results of the MC simulation 
and the correction factors /pb obtained from simulation, along with the estimated systematic 
uncertainties of the simulated values, are shown in Fig. [TJ The relative systematic uncertainties 
on P^e and P^ e b obtained by simulation have been estimated to be 10%, and are mainly due to the 
simulation of cluster reconstruction and energy calibration. However the error of the ratio /pb = 
Pfj-e/Pjie 1S significantly smaller (<5/pb//pb = 2%) due to cancellation of the main systematic 
effects. The measured Pj^ is in agreement with the simulation within their uncertainties. 

The K^2 background contamination integrated over lepton momentum has been computed 
to be (6.11 ± 0.22)% using the measured P^ b corrected by /pb- The quoted error is due to the 
limited size of the data sample used to measure Pj^ (0.16%), the uncertainty 5/pb (0.12%), and 
the model-dependence of the correction for the M^ iss (e) vs E/p correlation (0.08%). 

Rk is defined to be fully inclusive of intern al bremsstrahlung (IB) radiation^. The structure- 
dependent (SD) K + —> e + ^7 process D-7 | 18| ma y \ ea( \ to a j( e2 signature if the positron is 
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Figure 2: Left: Reconstructed squared missing mass M miss (e) distribution of the K e 2 candidates compared with 
the sum of normalised estimated signal and background components. Right: Lepton momentum distributions of 
the K e 2 candidates and the dominant backgrounds; the backgrounds are scaled for visibility. 



Table 1: Summary of backgrounds in the K e 2 sample. 



Source 


N B /N(K e2 ) 


K&in e ) 
K+ -► e+z/7 (SD+) 
K + -)■ ir°e + i> 
K+ -> vr+vr 
Beam halo 


(6.11 ±0.22)% 
(0.27 ± 0.04)% 
(1.07 ±0.05)% 
(0.05 ± 0.03)% 
(0.05 ± 0.03)% 
(1.16 ±0.06)% 


Total 


(8.71 ± 0.24)% 



energetic and the photon is undetected. In particular, the SD + component with positive photon 
helicity peaks at high positron momentum in the K + rest frame {E* ~ Mk/2) and has a 
similar branching ratio to K e 2- The background due to K + — > e + v^f (SD~) decay with negative 
photon helicity peaking at E* ~ M^/A and the interference between the IB and SD processes 
are negligible. The SD + background contribution has been estimated by MC simulation as 
(1.07±0.05)%, using a recent measurement of the K + — > e + v^ (SD + ) differential decay rate^l 
The quoted uncertainty is due to the limited precision on the form factors and decay rate, and 
is therefore correlated between lepton momentum bins. 

The beam halo background in the K e 2 sample induced by halo muons (undergoing \i — ^ e 
decay in flight or mis-identified) is measured using a data-driven method, by reconstructing 
K^2 candidates from a control K~ data sample collected with the K + beam dumped, to be 
(1.16 ±0.06)%. Background rate and kinematical distribution are qualitatively reproduced by 
a halo simulation. The uncertainty is due to the limited size of the control sample. The beam 
halo is the only significant background source in the K^2 sample. Its contribution is mainly at 
low muon momentum, and has been measured to be (0.38 ± 0.01)% using the same technique 
as for the K e 2 sample. 

The numbers of selected K e 2 and candidates are 59813 and 1.803 x 10 7 , respectively (the 
latter samples has been pre-scaled by a factor of 150 at the trigger level). Backgrounds in the 
K e 2 sample integrated over lepton momentum are summarised in Table [TJ the total background 



Table 2: Summary of the uncertainties on Rk- 



Source 


5R K x 10 5 


Statistical 


0.011 


K^2 background 


0.005 


K + — > e + z^7 (SD + ) background 


0.001 


K + — > 7r°e + i/, K + —7- 7r + 7r° backgrounds 


0.001 


Beam halo background 


0.001 


Helium purity 


0.003 


Acceptance correction 


0.002 


Spectrometer alignment 


0.001 


Positron identification efficiency 


0.001 


1-track trigger efficiency 


0.002 


LKr readout inefficiency 


0.001 


Total systematic 


0.007 


Total 


0.013 



contamination is (8.71 ± 0.24)%, and its uncertainty is smaller than the relative statistical 
uncertainty of 0.43%. The M^ iss (e) and lepton momentum distributions of K e 2 candidates and 
backgrounds are shown in Fig. [2j 

The ratio of geometric acceptances A(Ku2) I 'A{K e 2) in each lepton momentum bin has been 
evaluated with MC simulation. The radiative K + — > e + u^f (IB) process, which is responsible for 
the loss of about 5% of the K e 2 acceptance by increasing the recon struct ed M^ iss (e), is taken 
into account following^^, with higher order corrections according to ^^ l 

The acceptance correction is strongly influenced by bremsstrahlung suffered by the positron 
in the material upstream of the spectrometer magnet (Kevlar window, helium, DCHs). This 
results in an almost momentum- independent loss of K e 2 acceptance of about 6%, mainly by 
increasing the reconstructed M^ iss (e). The relevant material thickness has been measured by 
studying the spectra and rates of bremsstrahlung photons produced by low intensity 25 GeV/c 
and 40 GeV/c electron and positron beams steered into the DCH acceptance, using special data 
samples collected in the same setup by the NA48/2 experiment in 2004 and 2006. The material 
thickness during the 2007 run has been estimated to be (1.56 ± 0.03)%A"o, where the quoted 
uncertainty is dominated by the limited knowledge of helium purity in the spectrometer tank. 

A x 2 ffi to the measurements of Rk in the 10 lepton momentum bins has been performed, 
taking into account the bin-to-bin correlations between the systematic errors. The uncertainties 
of the combined result are summarized in Table [2J To validate the assigned systematic uncer- 
tainties, extensive stability checks have been performed in bins of kinematic variables and by 
varying selection criteria and analysis procedures. The fit result is^H 

R K = (2.487 ± 0.011 stat . ± 0.007 syst .) x 10~ 5 = (2.487 ± 0.013) x 10~ 5 , (5) 

with x 2 /ndf = 3.6/9. The individual measurements with their statistical and total uncertainties, 
the combined NA62 result, and the new world average are presented in Fig. [3l 

4 Search for lepton number violation 

The — > 7r T /x =1: // ± decay has been searched for using the NA48/2 2003-04 data sample, nor- 
malizing to the abundant — > ■k^it^it^ normalization channel (denoted below). Three- 
track vertices (compatible with either — > -Kfi^L or decay topology) are reconstructed by 
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Figure 3: Left: Measurements of Rk in lepton momentum bins with their uncorrelated statistical uncertainties 
and the partially correlated total uncertainties; the average Rk and its total uncertainty are indicated by a band. 
Right: The new world average including the present result. 



extrapolation of track segments from the spectrometer upstream into the decay volume, tak- 
ing into account the measured Earth's magnetic field, stray fields due to magnetization of the 
vacuum tank, and multiple scattering. The vertex is required to have no significant missing 
momentum, and to be composed of one ir^ candidate (with the ratio of energy deposition in the 
LKr calorimeter to momentum measured by the spectrometer E/p < 0.95, which suppresses elec- 
trons, and no in-time associated hits in the MUV), and a pair of ^ candidates (with E/p < 0.2 
and associated signal in the MUV). The muon identification efficiency has been measured to be 
above 98% for p > 10 GeV/c, and above 99% for p > 15 GeV/c. 

The invariant mass spectra of the reconstructed tt^^^ and ir^^^ candidates are pre- 
sented in Fig. |U The observed flavour changing neutral current — > tt^/j.^^ decay (3120 
candidates with a background of 3.3%) has been studied separately^! In the mass spectrum 
with same sign muons, corresponding to the lepton number violating signature, 52 events are 
observed in the signal region \M n ^ - M K \ < 8 MeV/c 2 . The background comes from the 
decay, and has been estimated by MC simulation to be (52.6 ± 19.8) events. The quoted 
uncertainty is systematic due to the limited precision of MC description of the high-mass re- 
gion, and has been estimated from the level of data/MC agreement in the control mass region 
of (465; 485) MeV/c 2 . This background estimate has been cross-checked by fitting the mass 
spectrum in the region between 460 and 520 MeV/c 2 , excluding the signal region between 485 
and 502 MeV/c 2 , with an empirical function similar to that used in the E865 analysis^! using 
the maximum likelihood estimator and assuming a Poisson probability density in each mass bin. 

The Feldman-Cousins method^S is employed for confidence interval evaluation; the system- 
atic uncertainty of the background estimate is taken into account. Conservatively assuming 
the expected background to be 52.6 — 19.8 = 32.8 events to take into account its uncertainty, 
this translates into an upper limit of 32.2 signal events at 90% CL. The geometrical accep- 
tance is conservatively assumed to be the smallest of those averaged over the — >• tt^u^jj,^ 
and samples (A^^^ = 15.4% and A^ = 22.2%). This leads to an upper limit of 
BR(^T ± -> ^ T /i ± / u ± ) < 1.1 x 10" 9 at 90% CL, which improves the best previous limit ^ 
by almost a factor of 3. 
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Figure 4: Reconstructed M,r W spectra for candidates with different (left) and same sign (right) muons: data 
(dots), Kzt: and K n/ltl MC simulations (filled areas); fit to background using the empirical parameterization as 
explained in the text (solid line). The signal region is indicated with arrows. 



5 Conclusions 

The most precise measurement of lepton flavour violation parameter Rx has been performed: 
Rk = (2.487 ±0.013) x 10 -5 is consistent with the SM expectation, and can be used to constrain 
multi-Higgs^l and fourth generation^ new physics scenarios. An improved upper limit of 1.1 x 
10~ 9 for the branching fraction of the lepton number violating — > 7r =F // ± // = ' = decay has been 
established. 
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